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ANTENNA  OPTIMIZATION  CRITERIA 


David  K.  Qicng 


ABSTRACT 

This  is  the  final  report  for  Contract  No.  F30602-68-C-0067  (ARPA 
Order  No.  1010)  monitored  by  the  Rome  Air  Development  Center  The  effec¬ 
tive  period  of  this  Contract  was  from  22  September  1969  to  22  March  1971. 

lour  technical  reports  have  been  Issued  on  the  results  obtained 
under  this  Contract.  They  are: 

(1)  Tech.  Rpt.  No.  1:  "Spacing  Perturbation  Techniques  for  Array 
Optimization,"  RADC-TR-68-19,  November  1967. 

(2)  Tech.  Rpt.  No.  2;  "Array  Optimization  Criteria,"  RADC-TR-68-579, 
November  1968. 

(3)  Tech.  Rpt.  No.  3:  "Beam  Synthesis  Techniques  for  Large  Circular 
Arrays  with  Many  Directive  Elements,"  RADC-TK-69-4 11 ,  October  1969. 

(4)  Tech.  Rpt.  No.  4:  "Sidelobe-Reduction  and  Interference-Suppression 
Techniques  for  Phased  Arrays  Using  Digital  Phase  Shifters," 

RADC-TK- 70-51,  February  1970. 

The  present  report  is  divided  into  two  parts.  Part  (A)  summaries  tlie 
essential  optimization  techniques  for  antenna  arrays,  including  a  method 
which  makes  only  phase  adjustments.  Part  (B)  presents  a  new  integral - 
equation  approach  for  optimizing  arrays  with  mutual  coupling.  This 
.approach  is  particularly  advantageous  when  an  array  contains  many  long 
dipole  elements. 
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PART  (A).  OPTIMIZATION  TECHNIQUES  FOR  ANTENNA  ARRAYS 


I.  INTRODUCTION 


An  antenna  is  an  essential  part  of  any  electronic  system  wlilch  transmits 
or  receives  electromagnetic  energy  in  a  wireless  fashion.  Without  an  antenna, 
electromagnetic  energy  will  be  localized,  and  interaction  at  a  distance  be¬ 
tween  unconnected  points  in  space  does  not  occur.  An  antenna  can  be  considered 
as  a  transducer  which  converts  electromagnetic  waves  in  space  to  current  or 
voltage  variations  in  a  circuit,  or  vice  versa.  It  must  be  an  efficient  radi¬ 
ator  (or  collector)  of  electromagnetic  energy,  and  it  should  direct  the  energy 
to  certain  desired  directions  and  suppress  it  in  other  specified  directions. 

Thus,,  one  must  be  concerned  not  only  with  the  conversion  efficiency  of  an  an¬ 
tenna,  but  also  with  its  spatial  response,  or  radiation  pattern.  In  an  en¬ 
vironment  which  does  not  Involve  nonlinear  media,  a  reciprocity  relation  holds 
such  that  the  properties  (pattern,  gain.  Impedance)  of  an  antenna  used  for  re¬ 
ceiving  are  identical  with  those  when  it  is  transmitting.  For  simplicity  we 
shall  then  refer  all  discussions  to  radiation  properties- 

A  radiating  element  of  electromagnetic  energy  may  take  many  different 
forms.  It  may  be  a  piece  of  conducting  wire,  a  dielectric  rod,  a  metallic 
horn,  or  a  slot  on  t^e  side  of  a  waveguide.  The  radiation  pattern  of  a 
single  element  is  fixed  for  a  given  frequency  of  excitation  and  contains,  in 
general,  a  main  beam  and  a  number  of  smaller  sidelobes.  In  practical  appli¬ 
cations  there  is  quite  often  a  need  for  either  improving  the  directive  properties 
or  controlling  the  sldelobe  structure  of  the  radiation  pattern.  Two  methods  are 
available  for  this  purpose:  one  method  is  to  use  an  appropriately  shaped  re¬ 
flector  or  lens  fed  by  a  radiating  element,  and  the  other  is  to  employ  a 
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number  of  radiating  elements  properly  arranged  In  space  to  form  an  antenna 
array.  When  It  Is  necessary  to  steer  (scan)  the  main  beam  of  the  radiation 
pattern,  the  requisite  motion  of  heavy  reflector-or  lens>type  antennas  en¬ 
tails  both  mechanical  and  structural  problems.  Moreover,  the  possible  rate 
of  scan  is  severely  limited.  On  the  other  hand,  beam-steering  for  an  antenna 
array  can  be  accomplished  electronically  by  adjusting  the  relative  phase  of 
excitation  in  the  array  elements  with  no  need  for  mechanical  motion,  result¬ 
ing  in  a  phased  array.  We  shall  concern  ourselves  only  with  phased  array 
antennas  in  this  report. 

The  radiation  pattern  of  an  .<iirray  antenna  obviously  depends  on  both  the 
array  geometry  and  the  pattern  of  the  individual  array  elements  Aside  from 
such  circuit  properties  as  impedance  and  efficiency,  the  parameters  which 
characterize  antenna  performance  are  all  based  on  the  shape  of  the  radiation 
pattern.  Performance  optimization  then  is  a  procedure  for  the  maximization  or 
minimization  of  certain  measures  on  the  radiation  pattern.  One  well-known  re¬ 
sult  in  this  regard  is  the  Chebyshev  array  which  makes  use  of  the  properties 
of  Chebyshev-Akhiezer  polynomials  [1],  (2).  A  Chebyshev  array  is  optimum  in 
the  sense  that  for  a  specified  sidelobe  level  the  width  of  the  main  beam  of 
its  radiation  pattern  is  a  minimum.  Conversely,  for  a  specified  beamwidth 
all  Che  sidelobes  of  a  Chebyshev  array  are  of  equal  height  and  are  at  a  lowest 
level  The  original  Chebyshev  design  considered  by  Dolph  (1]  was  limited  Co 
linear  broadside  arrays  of  isotropic  elements  uniformly  spaced  at  a  distance 
equal  to  or  larger  Chan  a  half-wavelength.  It  has  since  been  extended  in 
various  ways  (see,  for  instance,  [3]).  Recently  methods  for  determining  Che 
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current  distribution,  the  minimum  required  number  of  controlled  elements, 
and  other  properties  of  optimum  rectangular  arrays  with  a  steerable  main 
bean  and  constant  sldelobes  have  been  formulated  [4,17]. 

Besides  the  beanwldth-sldelobe  relationship,  an  Important  performance 
Index  for  any  antenna  Is  Its  directive  gain,  or  directivity  Directivity  is 
defined  as  the  ratio  of  the  radiation  Intensity  (radiated  power  per  unit 
solid  angle)  In  the  direction  of  the  main  beam  to  the  average  radiation 
Intensity.  To  put  It  another  way,  the  directivity  of  an  antenna  or  an  array 
Is  the  ratio  of  Its  maximum  radiation  Intensity  to  the  radiation  intensity 

A 

of  an  Isotropic  (omnidirectional)  source  radiating  the  same  total  power. 

It  measures  the  ability  of  concentrating  the  radiated  energy  In  the  main- 
beam  direction.  Our  attention  In  this  report  will  be  directed  toward  the 
various  techniques  for  maximizing  the  directivity  of  antenna  arrays. 

The  performance  of  an  antenna  system  as  a  receiving  device  Is  often 
constrained  by  the  presence  of  a  spatially  distributed  background  noise  as 
well  as  by  the  noise  generated  in  the  receiving  system.  A  useful  perfor¬ 
mance  index  of  a  receiving  array  is  the  slgnal-to-nolse  power  ratio  (SNR) 
at  the  system  output.  The  problem  of  finding  the  complex  weighting  fac¬ 
tors  of  the  Individual  array  elements  such  that  the  SNR  is  maximized  for  a 
signal  coming  from  a  given  direction  and  a  noise  of  a  given  power  spectral 
density  and  a  given  spatial  distribution  Is  of  considerable  importance.  It 

A 

It  should  be  noted  that  an  Isotropic  source  radiating  uniformly  in  all 
directions  Is  physically  unrealizable  for  vector  fields. 
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can  be  shown  [5,6]  that  techniques  similar  to  those  for  the  maximization  of 
directivity  are  applicable  for  SNR  maximization.  The  formulation  becomes 
more  involved  when  the  element  space-frequency  response,  the  signal  power 
spectral  density,  the  internal  noise  pov/er  spectral  density,  the  spatial  noise 
cross-power  spectral  density,  and  the  receiver  or  filter  frequency  response  are 
to  be  considered  [7].  We  shall  not  attempt  to  discuss  the  more  general  case  in 
this  report. 

The  problem  of  maximizing  the  directivity  of  a  linear  array  with  equally 
spaced  isotropic  elements  was  first  studied  by  Uzkov  in  1946  [8].  He  demon¬ 
strated  that  the  maximum  obtainable  directivity  for  an  array  with  N  elements 

2 

spaced  at  a  half-wavelength  (A/2)  apart  is  N  and  that  it  tends  to  N  as  the 
spacings  approach  zero.  Bloch,  Medhurst  and  Pool  [9]  examined  the  maximum 
directivity  of  a  linear  array  of  half-wave  dipoles  from  the  point  of  view  of 
self  and  mutual  resistances  of  the  elements.  Gain  optimization  under  a  speci¬ 
fied  constraint  was  Investigated  by  Uzsoky  and  Solymar  [10],  and  Lo,  Lee  and 
Lee  [5].  In  1964  Tal  [11]  published  many  interesting  curves  showing  the  opti¬ 
mum  directivity  of  various  types  of  uniformly  spaced  broadside  arrays,  linear 
arrays  with  maximum  radiation  in  the  direction  of  the  array  normal.  The 
optimization  problem  was  generalized  by  Cheng  and  Tseng  [12],  [13]  to  Include 
arrays  of  non-isotroplc  elements  arranged  in  an  arbitrary  configuration  with 
a  main  beam  pointing  at  an  arbitrary  direction.  By  making  use  of  a  theorem  on 
the  properties  of  a  ratio  of  two  Hermitlan  forms  in  matrix  algebra,  the  optimi¬ 
zation  procs?dure  was  formalized  in  a  concise  manner-  It  turned  out  that 
Krupltskli  in  U.S.S.R.  [14]  had  used  the  same  theorem  to  prove  the  existence 


-4- 


and  uniqueness  of  a  solution  for  exciting  an  array  of  discrete  radiators  for 
maximum  directivity. 

In  the  following  we  shall  first  express,  In  Section  11,  the  directivity 
and  slgnal-to-nolse  power  ratio  of  an  array  of  discrete  elements  as  a  ratio 
of  two  Hermltlan  forms.  The  optimization  principle  for  a  ratio  of  Hermltlan 
forms  Is  then  reviewed  and  applied  to  antenna  arrays  In  Section  111.  With  a 
given  array  configuration  where  the  element  positions  are  not  to  be  changed, 
the  excitation  amplitudes  and  phases  In  the  array  elements  can  be  adjusted 
for  the  optimization  of  a  performance  Index.  If  the  array  has  a  total  of  N 
elements,  the  optimization  procedure  Involves  the  determination  of  2N  parame¬ 
ters.  Typical  results  for  linear  and  circular  arrays  will  be  presented.  For 
linear  arrays  the  spaclngs  between  the  array  elements  represent  another  con¬ 
venient  set  of  parameters  that  can  be  adjusted  to  Improve  the  performance 
Index  further.  A  spacing-perturbation  technique  which  can  be  used  In  con¬ 
junction  with  the  adjustments  In  excitation  amplitudes  and  phases  Is  dis¬ 
cussed  In  Section  IV.  This  process  of  excitation  and  spacing  adjustments 
may  be  repeated  until  the  Improvement  obtained  by  further  adjustments  is  no 
longer  significant. 

In  practice  It  Is  perhaps  Inconvenient  to  adjust  the  element  spaclngs. 
Even  adjustments  In  excitation  amplitudes  are  difficult  and  expensive  to 
make.  Techniques  for  maximizing  the  directivity  of  a  fixed  array  by  keeping 
the  amplitudes  equal  and  adjusting  the  phases  only  are  therefore  of  Interest. 
These  techniques  are  explored  In  Section  V. 
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The  adjustment  of  the  excitation  amplitudes,  phases,  or  spaclngs  of 
an  array  in  order  to  achieve  a  maximum  directivity  changes  the  array  radi¬ 
ation  pattern.  In  practice  it  is  often  desirable  to  control  some  aspects 
of  the  array  pattern.  For  example,  one  may  wish  to  have  a  maximum  direc¬ 
tivity  in  one  direction  while  requiring  a  null  in  certain  other  directions 
in  order  to  minimize  Interference.  The  method  of  optimization  under  con¬ 
straints  is  discussed  in  Section  VI. 

Initially,  as  we  develop  the  optimization  procedure,  we  shall  assume 
that  all  the  elements  in  an  array  are  identically  polarized  and  have  the 
same  radiation  pattern  (element  pattern).  Once  the  element  patteni  is 
specified,  the  physical  structure  of  the  array  elements  is  no  longer 
important  in  our  problem,  and  it  is  immaterial  whether  the  elements  are 
dipoles,  horns,  slots,  or  other  apertures.  This  assumption  neglects  the 
implications  of  mutual  coupling.  For  large  arrays  with  many  elements  this 
assumption  gives  acceptable  results,  although  the  element  pattern  must  first 
be  found.  However,  for  arrays  with  a  small  number  of  closely  spaced  elements 
or  in  cases  where  more  accurate  results  are  desired,  mutual-coupling  effects 
must  be  taken  into  account.  Section  VII  reviews  the  moment  method  for  opti¬ 
mizing  the  directivity  of  arrays  of  wire  antennas  without  neglecting  mutual 
coupling.  The  moment  method  provides  numerical  solutions  by  first  converting 
the  governing  Integro-dlfferentlal  equations  into  matrix  equations- 

Finally,  in  Section  VIII,  we  discuss  various  other  factors  which  are 
relevant  in  the  optimization  of  discrete  antenna  arrays. 
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II.  GENERAL  FORMULATION 


Consider  an  array  of  N  discrete,  similarly  oriented,  identical  elements 

arranged  arbitrarily  In  a  3-dlmenslonal  space,  as  shown  in  Fig.  1.  Denoting 

the  excitation  In  the  nth  element  located  at  (r  ,6  ,41  )  by  I  exp(jij/  ),  we  may 

n  n  n  n  n 

write  the  array  factor  for  electric  field  Intensity  as 


where 


N 

E(e,4i)  “  I  exp(J(4i^  +  kr^  cos  a^)  ] 
n"l 


cos  a  «  sin  6  sin  6  cos(4>-4>  )  -f  cos  0  cos  6 

n  n  n  n 


(1) 

(2) 


k  ■  2^/A  Is  the  wavenumber,  and  A  Is  the  operating  wavelength.  Let  g(6,4) 
denote  the  element  power-pattern  function  which  Is  normalized  such  that 


g(eo**^o^  ■  ^ 

In  the  direction  of  the  main  beam.  The  directivity  of  the  array 

Is  then 


Radiation  Intensity  In  direction  (6  ,((>  ) 

00 

Average  radiation  Intensity 

_ _ 

^2n  IT 

^  I  d4»  [  |E(0,4i)j^  g(0,4i)  sin  Od') 

0  6 


(A) 


We  now  define  two  N-element  column  vectors  J  and  F^,  with  J  repre¬ 
senting  the  set  of  complex  element  excitation  functions: 
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1^  exp(jii/^) 

I2  exp(Ji;/2) 


exp(Jii/jj) 


(5) 


and  representing  Che  set  of  phase  factors  due  to  differences  in  distance: 


F 

o 


on 


exp(-jkrj  cos 
exp(-jkr2  cos  0^2) 


exp(-jkr^  cos  , 


(6) 


where  cos  a  (n*lt2 . N)  Is  obtained  from  (2)  by  setting  0  ■  6  and 

on  o 

ti  •  fi  .  From  (1),  (A),  (5)  and  (6)  it  Is  readily  verified  that  the  direc- 
o 

civlty  can  be  written  as 

j"  A  J 

(7) 

J  B  J 

where  t  on  a  matrix  indicates  the  adjoint,  or  Che  conjugate  transpose,  of 
the  matrix.  A  and  B  are  N  by  N  square  matrices  defined  as 


and 

with 

b 

mn 


A 


(a_nl 

mn 


F' 


~o 


mn 


2’r  TT 
I  r  f 

•7—  I  di{i  I  g(0,4))  exp(-jk(r  cos  a  -r  cos  a  )]  sin  OdO 
An  J  r  I  ex  ir/  J  '  m  n  n 

0  0 


(8) 

(9) 


(10) 
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It  Is  obvious  that  matrices  A  and  B  are  both  Hermltlar,  l.e.,  A'  »  A 
a  j-  a 

(a  ■  a  ),  and  B‘  ■  B  (b  ■  b  ).  Hence  D  In  (7)  is  a  ratio  of  two  Hcrniltlan 
nm  mn  ’  ~  ~  nm  on 

forms. ^  In  addition,  B  is  positive-definite,  which  implies  that  for  any  £  0, 

J^B  J  >  0.  This  is  proved  in  reference  [13].  Since  the  elements  of  matrices 
A  and  B  are  known  when  the  array  geometry,  the  operating  wavelength  and  the 
scan  angle  are  given,  the  optimization  problem  reduces  to  the  determination  of 
the  excitation  matrix  J  such  that  D  in  (7)  is  maximized. 

In  receiving  systems  the  output  slgnal-to-nolse  ratio.  Instead  of  the 
directivity,  is  of  interest.  The  output  SNR  may  be  defined  as  the  ratio  of  the 
power  received  per  unit  solid  angle  in  the  direction  of  the  signal  to  the  aver¬ 
age  noise  power  received  per  unit  solid  angle.  If;  is  only  necessary  to  replace 
the  power-pattern  function  g(6,(ti)  in  (A)  by  a  more  general  weighting  function 
w(6,iti)  which  includes  the  spatial  distribution  jf  noise  power.  We  write 


w(e,((i)  -  g(e,«(.)  T(e,^)  ,  (ii) 

where  T(6,((i)  is  the  spatial  distribution  function  of  noise  power.  Here  we 
understand  noise  to  be  a  combination  of  Interference,  clutter,  atmospherics, 
and  random  noise.  It  is  cle-ar  that  replacing  g(6,(ti)  by  w(6,4>)  does  not  change 
the  nature  of  the  cptlmlzatlon  problem.  In  fact,  the  expression  for  SNR  reduces 
to  that  for  D  when  T(6,(ti)  ■  1.  We  expect  a  suppression  of  the  sldelobes  in  the 
directions  of  high  noise  power  for  SNR  Improvement  [6].  In  the  next  section 
we  state  the  theorem  for  maximizing  D  by  adjusting  J. 


Sometimes  Hermltlan  forms  are  written  as  general  inner  products;  for  instance. 


N  N 

j  rA  J  -  y  y  j"  a  J 
.v^A.  t-  t-  mmnn 


m-1  n«l 


-J,  A 


They  are  quadratic  forms  in  Hilbert  space  in  the  variables  J^. 
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III.  OPTIMIZATION  BY  EXCITATION  ADJUSTMENTS 


A  theorem  In  matrix  algebra  on  the  properties  of  a  ratio  of  two  Hermltlan 
forms  [15]  Is  useful  for  the  maximization  of  directivity  by  excitation  adjustments. 
It  may  be  stated  as  follows: 

Theorem  1  -  If  a  quantity  D  Is  expressible  as  a  ratio  of  two  Hermltlan 
forms  as  In  (7)  and  If  B  Is  nonsingular  and  positive  definite,  then 
the  largest  eigenvalue  of  the  "regular  pencil"  of  matrices  A  XB 

Is  the  maximum  obtainable  D  when  J  Is  the  eigenvector  satisfying  the 
homogeneous  equation 


A  J  -  X„  B  J 

^  ^  M  ^ 


(12) 


For  our  case,  the  following  corollary,  proved  in  [13],  makes  the  optimization 
procedure  particularly  straightforward  and  simple. 

Corollary  -  If  A  in  (7)  is  expressible  in  the  form  of  (8),  then 

(a)  the  largest  and  only  nonzero  eigenvalue  of  the  regular 
pencil  A  -  XB  is 

•  *v  ‘V 

M  M  -vO  O 

and 

(b)  the  eigenvector  corresponding  to  X^  Is 


in 


b"^f 

~o 


(14) 


Equations  (13)  and  (14)  solve  the  problem  of  optimization  by  excitation 
adjustments  for  an  arbitrary  array.  For  a  given  array  configuration,  it 
Is  only  necessary  to  determine  the  elements  of  the  matrices  F  and  B,  In 
accordance  with  (6)  and  (10)  respectively. 
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For  a  linear  array  with  elements  located  arbitrarily  at  distances  d 
from  a  reference  point, 


n 


a  =  exp[-jk(d  -  d  )sln  e  ]  (15) 

mn  n  m  o 

where  6  denotes  the  main-beam  direction  measured  from  the  normal  to  the 
o 

array,  and 


b 

mn 


2tt  7t 

d4i  [  g(6,<(>)  expl-jk(d^-dj^)sln  e]de  , 
0  0 


(16) 


If  the  array  elements  are  Isotropic,  g(d,()>)  =  1,  and  are  equally  spaced, 

d  -d  “  (n-m)d,  (15)  and  (16)  reduce  to 
n  m 


and 


a  =  exp[1k(m-n)d  sin  6  ] 
mn  o 

_  sin  k(m-n)d 
mn  (m-n)d 


(17) 

(18) 


As  an  example,  It  has  been  shown  [12]  that  an  endflre  array  with  8 

Isotropic  elements  equally  spaced  at  0.425A  apart  has  a  directivity  of  12.5 

* 

with  a  uniform  amplitude  and  cophasal  excitation.  Optimization  by  the  above 
procedure  results  In  a  directivity  of  22.0.  The  amplitude  Is  tapered  (center- 
to-edge  ratio:  1.69)  and  the  phase  shift  between  adjacent  elements  Is  roughly 
170®.  Other  examples  of  linear-array  optimization  by  excitation  adjustments 


An  excitation  Is  said  to  be  cophasal  when  the  relative  phases  In  the  elements 
are  adjusted  such  that  the  contributions  of  all  the  elements  add  In  phase  In 
the  main-beam  direction.  Thus,  for  an  endflre  array  with  0.425A  spacing,  the 
progressive  phases  In  the  neighboring  elements  differ  by  0-425  >■  360*or  153°, 
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will  be  given  In  Section  IV  when  spacing-perturbation  techniques  are  dis¬ 
cussed. 

Besides  the  linear  array,  the  circular  array  represents  another  class 
of  arrays  with  a  simple  geometry  and  Important  practical  applications.  Figure 
2  shows  a  general  circular  array  with  nonunlformly  spaced  elements  In  the  xy-plane. 
Substituting  0^  =  ti/2  In  (3),  we  obtain 


cos  a  =  sin  6  cos  )  . 

n  n 


(19) 


Since  r  =  p  =  p  for  all  n,  b  In  (10)  simplifies  to 
n  n  '  mn 


2tt  tt 


b  =  ■— 
mn  4ti 


d{|i 


g(9,(J))  exp{-Jkp  cos(ij)-ij)  )  sin  6}  sin  0d6  , 

r  J  UJJJ  T  jjjjj 


(20) 


where 


and 


p  -  2p  sln((J»  -4  )/2| 
mn  m  n 


,  sin  i)>  -  sin  ij) 

^ _  *  tan  . . )  . 

mn  cos  ij)  -  cos 


(21) 


(22) 


m  n 

For  a  circular  array  of  N  uniformly  spaced  Isotropic  elements, 


g(0,({i)  =*  1,  ®  2nTi/N,  and 


p  *  2p  sin  (m-n)ii/N  I  . 
mn  ' 


(23) 


We  have,  from  (20), 


sin  (kp  ) 
mn 


mn  kp 


(24) 


mn 


The  expressions  of  b^^  for  some  directive  elements  with  typical  power- 


pattern  functions  have  also  been  given  [13]. 
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Figure  3  compares  the  maximum  directivity  with  the  directivity 
under  a  uniform-amplitude  and  cophasal  excitation  for  a  12-element 
circular  array  as  a  function  of  the  array  diameter.  We  note  that  is 
everywhere  higher  than  and  that  Increases  very  rapidly  when  the 
array  diameter  is  less  than  2A  (a  superdlrectlve  situation).  Since  super¬ 
directive  arrays  require  very  large  currents  of  opposite  signs  in  neighbor¬ 
ing  elements,  resulting  in  excessive  heat  loss  and  very  low  radiation  in¬ 
tensity  in  the  direction  of  the  main  beam,  it  is  appropriate  to  define  a 
main-beam  radiation  efficiency  n.  A  practical  optimum  design  then  would 
be  a  suitable  compromise  between  and  n.  We  define 


|E(9  ^  )| 

n  - - 2 - ° - -  100%  . 


N 


(25) 


Nil, 

n-1 


n 


By  the  use  of  Schwarz's  inequality,  it  is  easy  to  show  that  [13]  n  equals 
100%  only  for  uniformly  excited  cophasal  arrays  and  becomes  very  small  under 
superdlrectlve  situations.  The  main-beam  radiation  efficiency  turns  out  to  be 
the  reciprocal  of  the  tolerance  sensitivity  used  by  Uzsoky  and  Solymar  [10]  to 
measure  the  mean-square  variation  of  the  maximum  field  with  respect  to  the 
mean-square  deviation  of  the  excitation.  The  values  of  n  under  the  condi¬ 
tions  for  are  also  plotted  in  Fig.  3. 

It  has  been  pointed  out  [16]  that  the  optimization  problem  becomes 
easier  by  applying  an  orthogonalizatlon  process  for  arrays  possessing  a  cyclic 
symmetry.  For  large  circular  arrays  with  many  uniformly  spaced  elements,  the 
tedium  of  inverting  a  large  B  matrix  (in  (14))  can  be  circumvented  through  the 


introduction  of  a  rotational  operator  [31]. 
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IV.  OPTIMIZATION  BY  SPACING  ADJUSTMENTS 


In  Che  preceding  section  optimization  was  achieved  by  adjusting  the 

excitation  amplitudes  and  phases  In  the  fixed  elements  of  a  given  array.  If 

Che  element  positions  are  also  allowed  to  vary,  we  acquire  an  additional 

dimension  of  freedom  (which  represents  an  additional  N-1  degrees  of  freedom 

for  a  linear  array  with  N  elements),  and  we  expect  to  be  able  to  Improve  on 

Che  results  obtained  by  excitation  adjustments  only.  A  spacing-perturbation 

technique  exists  which  Is  useful  for  optimizing  the  directivity  In  a  given 

direction  or  Che  slgnal-to-nolse  ratio  In  a  given  noise  environment  [6]. 

This  technique  will  be  developed  In  this  section. 

Consider  a  linear  array  of  N  Identical  elements  symmetrically  located 

about  the  origin  along  the  x-axls.  Let  (6p»  be  the  direction  of  the 

main  beam,  and  I  exp(j4>  )  and  I  exp(-j4>  )  be  the  excitations  In  the  nth  and 
n  n  n  n 

-nth  elements  respectively,  where 

$  =  kd  sin  9  cos  <P  +  \p  .  (26) 

n  n  0  on 

In  (26),  d^  Is  the  distance  of  the  nth  element  from  the  origin,  and  Is 
the  phase  shift  from  cophasal  operation.  Note  that,  except  for  the  assumed 
symmetry  about  the  origin,  there  Is  no  restriction  on  element  spaclngs  which 
may  be  nonuniform  and  the  elements  themselves  need  not  be  omnidirectional. 

■k 

If  N  Is  odd  and  equals  2M-1,  the  array  factor  Is 

•k 

Only  very  minor  modifications  are  needed  when  N  Is  even.  The  formulation 
is  entirely  parallel. 
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(27) 


M 

E(u)  =  I  +2  T  I  co8(6  u  +  ) 

o  n  n  n 

n=l 

where 

u  =  kd(sln  6  cos  4i  -  sin  6^  cos  (Ji^)  (28) 


6 

n 


d  /d 
n 


(29) 


and  d,  a  normalizing  distance,  may  be  any  choice  of  convenience.  For  ex¬ 
ample,  If  one  starts  with  an  equally  spaced  array.  It  would  be  natural  to 
make  d  the  uniform  spacing  between  neighboring  elements.  We  write  the 


output  slgnal-to-nolse  ratio  In  the  direction  of  the  signal  (6  ■  6^,  (})  =  (fi^, 
u  a  0)  as 


(30) 


where  w(e,({i)  is  a  weighting  function  defined  in  (11).  Our  problem  is  to 
find  the  set  of  normalized  element  positions  {6^}  such  that  SNR  Is  maxi¬ 
mized  for  a  given  set  of  excitation  amplitudes  and  phases.  Note  that  this 
becomes  a  directivity  maximization  problem  when  w(e,ij))  =  g(e,4>). 

Let  {6°}  denote  the  original  normalized  element  positions,  and 


+  X 

n 


(31) 


where  x^  represents  the  spacing  perturbation  for  the  nth  element  and 
x^  <<  1.  Substitution  of  (31)  In  (27)  yields  approximately 

M 

E(u)  =  E°(u)  -  2  I  I  X  u  sln(6°u  +  \l>  )  (32) 

‘'i  n  n  n  n 

n=l 
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where  E°(u)  is  the  original,  unperturbed  array  factor  with  6°  substituted 
for  in  (27).  Using  (32),  we  can  write  (30)  in  the  following  form: 


SNR 


[E°(^ 


-  2x'6  +  X 


^ 


(33) 


where 


2-n  TT 

d<t) 


1e°(u)|^  w(d,^)  sin  6d0 


(34) 


x'  -  (x^,X2, 


*n» 


(35) 


is  the  transpose  of  the  column  matrix  of  spacing  perturbations  x;  ^  is  a 
column  matrix  of  typical  element 


2v  n 


Bn  “  ^  j  d(J)  j  I^uE°(u)w(6,4i) 


sin  (6  u  +  )  sin  6d6; 

n  ’ 


(36) 


and  C 


[C  ]  is  an  N^N  square  matrix  with 
an 


mn  TT 


2it  n 


d$ 


I  I  u  w(e,iji) 
ran 


sin  (6°u  +  ij;  )  sin  (6°u  +  ij;  )  sin  9de  ,  (37) 

ra  ra  n  n 

It  can  readily  be  shown  that  C  is  symmetric  and  positive  definite. 
Use  can  then  be  made  of  the  following  theorem  which  is  proved  in  the 
Appendix. 
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Theorem  2  -  If  a  quantity  SNR  can  be  expressed  in  terms  of  an  N^l 


real  colunn  vector  x  as  in  (33),  where  a  is  a  constant,  6  is  another 
Nxl  real  colunn  vector,  and  C  is  an  N^n  positive  definite,  symmetric, 
square  matrix,  then 

Max  SNR  -  ■  :  (38) 

£'  iS 

and 

”9) 

Equations  (38)  and  (39)  give  the  results  of  a  first-order  per¬ 
turbation.  After  the  components  of  x^  have  been  determined  from  (39), 
one  can  then  use  (6°  +  }^)  as  the  new  normalized  element-position  column 
matrix  and  perform  a  second-order  perturbation  to  obtain  further  improvement 
in  the  performance  index.  This  process  can  be  repeated  until  it  becomes 
evident  that  further  iteration  yields  a  negligible  improvement.  The  final 
values  of  (6^)  determine  the  element  positions  for  a  maximum  SNR  for  the 
given  excitation.  Now  this  perturbed  nonuniformly  spaced  array  can  be 
further  optimized  by  proper  amplifications  and  phase  shifts  in  the  array 
elements  using  the  method  developed  in  the  preceding  section  [6],  A 
second  local  maximum  will  be  reached,  which  may  possibly  be  further  im¬ 
proved  by  holding  the  excitation  unchanged  and  again  perturbing  the  spacings. 
The  cycle  may  be  repeated  until  further  adjustments  are  no  longer  worthwhile. 

We  illustrate  the  application  of  the  above  technique  with  a  broadside 
array  of  7  isotropic  elements.  It  is  desired  to  optimize  the  array  for 
maximum  SNR  in  a  noise  environment,  T(u)  specified  by  Fig.  A  with  a>15. 
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■  (1/4)  (2'’d/\)  and  U2  *  (1/12)  (2nd/A) .  The  normalizing  distance,  d, 

Is  chosen  Co  be  0.8851  which  corresponds  to  Che  spacing  for  maximum  dlrec- 
clvicy  In  a  uniformly  spaced  linear  array  with  7  isotropic  elements.  The 
results  for  the  (a)  space-perturbed,  (b)  excitation-adjusted  (by  amplifi¬ 
cation),  and  (c)  optimized  arrays  are  listed  In  Table  1.  Ue  note  that 
large  improvements  in  SNR  are  possible  by  optimization  through  either 


Table  1. 


SNR  Optimization  for  Seven-Element  Broadside  Array 


1  .2nd. 

4  '  X 


_1_  .2nd. 
12  '  A 


H 

B 

B 

^3 

IRIi 

BBSS 

BBS 

SNR 

Uniform  array 

HQ 

1.00 

1.00 

1.00 

1.77 

1.77 

1.77 

16.0 

Space-perturbed 

array 

same  as  above 

1.65 

1.76 

2.10 

19.8 

Exc. -adjusted 
perturbed  array 

1.00 

0.89 

0.67 

0.39 

same 

as  above 

78.1 

Optimized  array 

1.00 

0.86 

0.59 

0.40 

1.66 

_ 1 

1.72 

1.74 

181.9 

spacing  perturbation  or  excitation  adjustments,  and  that  the  SNR  of  the 
optimized  array  is  about  11.4  times  that  of  the  original  uniform  array. 

Of  course,  Che  array  optimized  for  a  maximum  SNR  is  not  the  same  as  the  one 
for  a  maximum  directivity.  The  directivity  of  the  SNR-optimlzed  array  in 
Table  1  (sketched  in  Fig.  5(a))  is  10.6,  whereas  a  directiviry-optlmlzed 
broadside  array  of  7  isotropic  elements  has  a  directivity  of  10.63  [6]. 
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The  radiation  patterns  of  the  SNR-optimized  array  are  plotted  in  Fig.  5(b), 
where  the  spatial  distribution  of  the  noise  or  interference  power  is  also 
shown.  It  is  Interesting  to  see  that,  at  the  expense  of  a  slightly  wider 
main  beam,  the  sidelobes  of  the  optimized  array  are  everywhere  lower  than 
those  of  the  uniform  array.  In  particular,  the  first  sidelobe,  which 
normally  occurs  in  a  region  where  the  noise  power  is  high,  is  much  sup¬ 
pressed  and  its  position  slightly  shifted. 
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V.  OPTIMIZATION  BY  PHASE  ADJUSTMENTS 


In  Section  III  we  discussed  the  method  for  maximizing  array  directivity 
by  adjusting  the  excitation  amplitudes  and  phases  In  the  array  elements.  Ac¬ 
curate  amplitude  adjustments  require  the  use  of  precision  amplifiers  or  attenu¬ 
ators,  or  specially  designed  directional  couplers.  It  Is  hence  of  both  theo¬ 
retical  and  practical  Interest  to  develop  a  technique  for  optimizing  the 
directivity  of  uniformly  excited  arrays  requiring  only  phase  adjustments.  The 
elimination  of  the  need  for  amplitude  adjustments  would  result  In  a  simplified 
feed  structure  and  a  reduced  cost. 

Directivity  maximization  under  the  constraint  of  a  tnlform  amplitude 
in  all  the  array  elements  can  be  formulated  by  Lagrange  multiplier  methods  In 
several  ways.  However,  It  was  found  that  the  resulting  equations  were  not 
amenable  to  a  stable  solution  even  by  Iterative  methods,  because  of  con¬ 
vergence  problems.  On  the  other  hand,  a  perturbation  procedure  similar  to 
that  employed  In  Section  IV  for  spacing  adjustments  can  be  used  for  the  phase- 
adjustment  problem  [18].  The  essential  steps  of  this  procedure  will  be 
developed  In  the  following. 

Consider  a  linear  arr:ay  of  2^H■1  symmetrically  located,  equally  spaced, 
identical  elements,  all  with  the  same  excitation  amplitude  I.  The  array 
factor  is  then,  from  (27), 

M 

E(u)  ■  I[1  +  2  I  cos(nu  +  )  J  ,  (40) 

n-1  " 

where  u  is  defined  in  (28)  and  Is  the  phase  shift  from  cophasal  operation. 

We  may  write 
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(41) 


{1/  "  li)  +  X 

n 


where  tli  Is  an  assumed  Initial  value  and  x  <<  1  Is  a  small  perturjatlon. 
n  n 

Ordinarily  it  is  convenient  to  start  from  a  cophasal  excitation,  i.e., 
Substituting  (41)  in  (40),  we  obtain  approximately 


M 

E(u)  •  E°(u)  -  21  y  X  sln(nu  +  ij(°)  . 

*■ ,  n  n 

n«l 

Using  (42)  in  (4),  we  can  express  the  directivity  of  the  array  in  the 
following  form: 


(42) 


a.  -  x’S,  +  x'C,  X 


(43) 


where 


X' 


[  ■  •••#  *  *2*  ***•  (44) 

and  a,,  6,,  and  C,  are  similar  to  a,  6,  and  C  in  (33).  a,  is  ldent:.cal 
to  a  in  (34),  and  the  elements  of  S,  and  C,  can  be  obtained  froiii  (36) 

"l  -vl 

and  (37)  respectively  with  I  ■  1  ■  I  and  the  argument  (6'^u  -f  i;^  )  of  the 

m  n  n  n 

sine  functions  replaced  by  (nu  (^°).  Element  power-pattern  function 
g(0,<|i)  replaces  weighting  function  w(9,iti)  in  computing  directivity. 

It  is  now  obvious  that  the  same  Theorem  2,  which  is  proved  in  the 
Appendix  and  found  useful  for  optimization  by  spacing  adjustments,  can  be 
used  for  maximizing  0  in  (43)  by  phase  adjustments.  The  required  phase 
changes  in  the  array  elements  are  determined  from 

-  cT^B,  (45) 

"'1  •>'1 

which  represents  a  first-order  perturbation  from  the  initial  values 

One  may  then  treat  {1;^°  +  x^}  as  the  new  initial  values  and  perform 
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a  second-order  perturbation.  This  process  may  be  repeated  until  It  Is 
apparent  that  a  maximum  directivity  has  been  obtained. 

Although  the  preceding  formulation  for  optimization  by  phase  adjust¬ 
ments  starts  with  a  linear  array.  It  Is  clear  that  the  procedure  can  be 
applied  to  an  arbitrary  three-dimensional  array.  In  particular,  for  a 
circular  array  of  N  elements  with  radius  p  In  the  xy-plane,  a  phase  per¬ 
turbation  as  Indicated  In  (41)  results  In  an  array  factor 

N 

E(6,4i)  =  E°(0,(|i)  -  21  y  X  sln(A  +  1/;°)  (46) 

n  n  n 

n*l 

with 


A  ■  kp[sln  0  cos((}i-(}»  )  -  sin  0  cos(i)>  -i)>  )]  (47) 

n  noon 

where  i^i  denotes  the  location  of  the  nth  element  and  (0  ,()>  )  Is  the  dlrec- 
n  0  0 

tlon  of  the  main  beam.  We  note  that  (46)  Is  entirely  similar  to  (42). 
Substitution  of  (46)  In  (4)  will  yield  a  directivity  expression  In  the  form 
of  (43),  and  hence  the  same  optimization  procedure  follows.  The  optimum 
directivity,  D^,  obtained  by  phase  adjustments  only  for  a  circular  array 
with  12  uniformly  spaced  short  dipoles  Is  plotted  In  Fig.  6  as  a  function 
of  array  diameter.  In  the  same  figure  are  also  plotted  the  maximum 
directivity  when  both  amplitudes  and  phases  are  adjusted,  and  D^,  the 
directivity  under  a  uniform-amplitude  and  cophasal  excitation.  We  note  that 
the  curve  lies  everywhere  between  the  and  curves.  For  most  array 
diameters  less  than  3a  (element  spacing  less  than  3A/4)  an  Improvement  of 
about  2  dB  In  directivity  Is  possible  by  phase  adjustments  alone.  When  the 
array  diameter  Is  very  small,  the  directivity  of  the  phase-adjusted  array 
increases  rapidly.  Indicating  a  superdlrectlve  situation  which  Is  absent  In  a 
uniform  cophasal  array.  The  main-beam  radiation  efficiency  of  a  superdlrectlve 
array  tends  to  be  very  low,  as  has  been  pointed  out  In  Section  III. 
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VI.  OPTIMIZATION  WITH  CONSTRAINTS 


In  previous  sections  we  discussed  techniques  for  maximizing  the 
directivity  or  the  slgnal-to-nolse  ratio  of  an  antenna  array  without  con¬ 
straints;  that  Is,  without  Imposing  at  the  same  time  a  requirement  on  any 
other  performance  index  of  the  array.  We  have  already  seen  that  a  maximum 
directivity  for  an  array  with  closely  spaced  elements  Is  accompanied  by  a 
low  main-beam  radiation  efficiency.  In  practice,  we  may  desire  that  a 
maximum  directivity  be  obtained  together  with  a  prescribed  value  for  main- 
beam  radiation  efficiency.  The  existence  of  constraints  cr  auxiliary  con¬ 
ditions  effectively  reduces  the  total  number  of  Independent  variables  which 
can  be  adjusted  for  optimization.  In  such  cases,  a  procedure  using  Lagrange 
multipliers  can  be  applied  to  determine  the  stationary  value  of  directivity 
(19]-  This  approach  has  been  employed  [5]  to  maximize  slgnal-to-nolse  ratio 
under  a  constraint  on  a  Q-factor.  The  Q-factor  Is  a  quantity  which  Is  pro¬ 
portional  to  the  ratio  of  the  directivity  and  the  main-beam  radiation  ef¬ 
ficiency.  It  turns  out  that  this  approach  results  In  a  rather  Involved 
numerical  procedure.  We  shall  not  go  Into  It  further  here. 

A  more  useful  class  of  optimization  problems  with  constraints  per¬ 
tains  to  the  maximization  of  some  performance  index  while  controlling  the 
array  pattern  In  certain  definite  ways.  For  Instance,  one  may  wish  to  have 
a  maximum  directivity  In  the  direction  of  some  distant  transmitter  or  receiver 
and,  at  the  same  time,  to  minimize  the  Interference  from  some  other  directions. 
This  Is  equivalent  to  the  problem  of  directivity  maximization  with  controlled 
locations  of  certain  pattern  nulls,  and  can  be  reduced  to  that  of  maximizing 
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the  ratio  of  two  Hermltlan  forms,  as  studied  In  Section  Ill,  through  the 

Introduction  of  a  constraint  matrix. 

The  array  factor  In  (1)  can  be  written  as  the  Inner  product  of  the 

space  vector  F  <■  [exp(-1kr  cos  a  )]  and  the  excitation  vector  J  defined 
/V  n  n  ‘v 

In  (5);  l.e., 

E(0.«)  -  <F,2>  --  F^J  .  (48) 

Pattern  nulls  In  the  directions  (6^,  ()>^)  are  specified  by  homogeneous 
equations 

+  ^ 

!l  ^  =  I  ^*^n  “in^^  “  ° 

n=l 

1  =  1,  2 . M  <  (N-1) 

where  are  obtained  from  (2)  by  setting  0  •  0^  and  ())  •  4)^.  A  geo¬ 
metrical  Interpretation  of  (49)  Is  that  the  excitation  vector  J  which 
we  seek  to  maximize  the  directivity  In  the  direction  (6^,  ij)^)  Is  now 
required  to  be  simultaneously  orthogonal  to  M  Independent  constraint 
vectors  F^.  The  N-dlmenslonal  space  Is  divided  Into  two  mutually  or¬ 
thogonal  subspaces:  an  M-dlmenslonal  subspace  containing  the  constraint 
vectors  and  an  (N-M) -dimensional  subspace  where  the  excitation  vector  J 
must  lie.  The  mathematical  procedure  [21]  for  maximization  under  con¬ 
straints  consists  of  (a)  finding  an  appropriate  set  of  M  mutually  or¬ 
thogonal  vectors  that  occupy  the  same  subspace  as  the  constraint  vectors 

F,,  Ft . F,,,  (b)  obtaining  an  additional  orthogonal  (N-M)  vectors 

which  at  the  same  time  are  orthogonal  to  the  first  M  vectors, (c)  forming 
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from  steps  (a)  and  (b)  an  N  b\  N  normalized  constraint  matrix  which  is 
the  unitary  matrix  for  coordinate  transformation,  (d)  transforming  the 
directivity  expression  as  a  ratio  of  two  Hermltlan  forms  to  the  new 
orthogonal  coordinate  system,  and  (e)  maximizing  the  directivity  In  the 
same  manner  as  outlined  In  Section  Ill. 

The  Gram-Schmldt  procedure  [19]  provides  a  method  for  determining 
an  orthonormal  basis  for  a  vector  space  In  which  any  set  of  spanning  vec¬ 
tors  Is  known.  This  method  can  be  used  to  find  the  constraint  matrix  In 
steps  (a)  and  (b) .  Assuming  P  to  be  the  normalized  constraint  matrix 
which  Is  also  the  transformation  matrix,  we  write 


P  J  -  J  . 


(50) 


The  directivity  In  (7)  becomes,  after  the  coordinate  transformation. 


j'*'  (P  A  P^)J  A  J 

m^Q  A#  ^  #vC  -vC  *^0 

(P  B  P^)J  B  J 

^  A*  C  C 


(51) 


Inasmuch  as  each  of  the  first  M  mutually  orthogonal  vectors  is  a  linear 

combination  of  the  constraint  vectors,  the  first  M  rows  of  J  In  (50)  are 

linear  cond)lnatlons  of  the  homogeneous  constraint  equations  (49)  and  are 

therefore  zeros.  Hence  the  first  M  entries  in  J  and  the  first  M  rows 

~c 

and  M  columns  of  A^  and  B^  can  be  discarded,  resulting  In  an  abridged  form 
for  directivity  [20]: 


J  A  J 
,  ~a  ~a  ~a 

^a  “  1 -  ’ 

J  B  J 
-a  ~a  -a 


(52) 


where  A  and  B  are  the  abridged  (N-M)  by  (N-M)  matrices  and  J  is  the 
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abridged  (N-M) -element  column  vector.  The  remainder  of  the  optimization 
process  then  follows  In  exactly  the  same  manner  as  that  pertaining  to  the 
unconstrained  D  In  (7),  Section  III,  except  that  the  numerical  problem  Is 
now  simpler  because  the  matrices  Involved  are  of  a  lower  rank  It  is 
obvious  that  the  maximum  obtainable  directivity  with  constraints  will  be 
less  than  that  without  constraints  on  account  of  the  reduced  freedom. 
Typical  results  on  maximum  directivity  with  null  placements  and  with  re¬ 
duced  radiation  level  In  an  angular  sector  have  been  published  [20]. 
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VII.  CONSIDERATION  OF  MUTUAL  COUPLING 


One  tacit  assumption  implied  in  the  optimization  techniques  which  we 
have  considered  thus  far  is  that  all  the  elements  in  an  array  have  the  same 
radiation  pattern.  This  is  tantamount  to  assuming  that  mutual-coupling 
effects  are  negllglbl'a.  For  large  arrays  with  many  elements  this  assump¬ 
tion  is  acceptable  if  our  interest  lies  in  array  directivity  and  not  in  the 
current  distribution  or  the  exact  radiation  pattern  of  each  element.  The 
consideration  of  mutual-coupling  effects  in  array  optimization  greatly  com¬ 
plicates  the  problem.  To  the  author's  knowledge,  no  work  has  been  published 
on  array  synthesis  or  optimization  for  mutually  coupled  aperture-type  radi¬ 
ators.  With  wire  antennas  the  method  of  moments  can  be  used  to  obtain 
numerical  answers  [23-26].  In  this  section  we  will  outline  the  optimization 
procedure  for  arrays  of  wire  antennas  when  mutual  coupling  is  not  neglected. 

The  moment  method  for  solving  electromagnetic  problems  consists  mainly 
of  three  steps;  namely,  the  formulation  of  the  governing  Integro-dlf ferentlal 
equations,  the  expansion  of  the  unknown  functions  in  terms  of  a  set  of  linearly 
Independent  basis  functions,  and  the  testing  of  the  expanded  equation  by  form¬ 
ing  inner  products  with  a  set  of  linearly  Independent  set  of  weighting  func¬ 
tions  [26].  The  result  is  a  set  of  simultaneous  equations  which  can  be 
solved  by  matrix  methods.  The  choice  of  the  basis  and  weighting  functions 
depends  on  the  desired  accuracy  and  the  ease  in  evaluating  the  coefficients 
in  the  simultaneous  equations.  One  convenient  technique  in  the  choice  of 
basis  functions  is  to  divide  the  domain  of  the  unknown  function  into  small 
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Intervals  or  subsections.  Simple  basis  functions  such  as  pulses  or  triangles 
are  defined  to  exist  over  one  or  a  few  such  subsections  and  to  be  zero  else¬ 
where.  The  simplest  weighting  functions  are  Dirac  delta  functions  defined 
at  discrete  points  at  which  the  expanded  governing  equations  are  to  be  satis-- 
fled.  This  Is  the  point-matching  or  sampling  method. 

For  thin  wire  antennas,  the  currents  and  charges  on  the  wires  can  be 
approximated  by  current  and  charge  filaments  along  the  antenna  axes.  We 
consider  an  array  of  thin  linear  antennas  parallel  to  the  z-axls.  From  the 
Maxwell's  equation  for  tlme-harmonlc  fields  In  a  homogeneous  medium 


7  «  H  -  jue  E  +  J 


and  H  ■  —  7  *  A,  we  have 


Now 


and 


E--~(-7x7xA-j). 
Jwe  'u 


7x7xA-77.A-7^A 


(7^  +  k^)  A  •  -  UJ, 


2  2 

where  k  ■  u  ue .  Combination  of  (55)  and  (56)  with  (5A)  yields 


(53) 

(54) 

(55) 

(56) 


E  -  7^  (7  7  •  A  +  k^A)  . 

jwue 


(57) 


Since  J  has  only  a  z-component,  we  can  rewrite  (57)  as  a  scalar  equation 

dZ 


_  1 _  .8“  ^  .  2v. 

Z  Ju)U£  •x-2  z 


(58) 


For  thin  linear  antennas,  the  tangential  electric  field  E  at  the  center 

zp 

of  a  typical  pth  subsection  Is  then  [25] 
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2p  jue  ^2 


all 

antennas 


I(z')exp[-Jk|?^  -  ?'| ] 


411  |?p  -  V 


dz'  , 


(59) 


where  Is  the  position  vector  to  the  center  of  the  pth  subsection  under 
consideration,  and  r'  Is  a  vector  from  the  origin  to  a  point  on  an  antenna 
at  which  the  current  Is  I(z').  If  each  of  the  N  antennas  In  the  array  Is 
divided  Into  S  subsections  of  length  which  carries  a  constant  current  Iq, 

(59)  can  be  approximated  by  a  summation: 


E 

zp 


M 


r  113^  2  f  exp[-Jk|r  -  r’ | ) 

y  I^>T^(^+k^)  - 2 -  dz' 

-1  3z^  . 


4it  I  r  - 


(60) 


In  (60),  M  ■  NS.  The  quantity  In  the  wavy  brackets  Is  the  electric  field 
at  the  center  of  subsection  p  due  to  a  unit  current  In  subsection  q  and 
can  be  written  as  Z  /Al,  where  equal  subsections  (A£  •  AP.)  are  assumed 

pq  q 

for  simplicity.  Thus, 


E 

zp 


(A£)  •  V 

P 


M 

I  z 

qil  P‘1 


(61) 


or.  In  matrix  form, 

V  -  Z  I  ,  (62) 

where  V  «  [V  ]  and  1  ■  [1  ]  are  M  by  1  column  matrices  and  Z  -  [Z  ] 

Is  an  M  by  M  square  matrix.  Z  may  be  called  a  generalized  Impedance 
matrix  and  depends  only  on  the  geometrical  configuration  of  the  array. 

In  the  terminology  of  the  method  of  moments,  (62)  Is  the  result  of 
(a)  using  pulse  basis  functions  each  of  which  exists  over  only  one  sub- 
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section,  (b)  using  an  integral-type  inner  product,  and  (c)  using  Dirac  delta 
functions  as  weighting  functions.  The  differencial  operator  in  (60)  nay  be 
approximated  by  a  second-order  difference  operator  and  computer  subroutines 
for  the  calculation  of  Z  are  available  [25].  In  a  radiation  problem  V  is 

pq 

a  column  matrix  of  known  voltages  which  are  all  zero  except  for  Che  exci¬ 
tation  voltages  at  feed  points.  Column  matrix  I  determines  the  current 
distributions  on  the  wire  antennas: 


I  •  Z'^V  -  Y  V  .  (63) 

With  £  known,  all  field  quantities  of  interest  can  be  determined. 

The  array  directivity  defined  in  (4)  can  be  expressed  in  terms  of  the 
total  power  input  Co  Che  array,  if  the  antenna  wires  are  assumed  to 

be  perfectly  conducting.  The  electric  field  in  the  far  zone  of  an  array 
consisting  of  z-directed  wire  antennas  is 


Thus , 


Eg  •  -  juAg  •  Ju»A^  sin  0 


D  - 


4iTr^|E(e^)|^/2 

2Tir^|a)A  sin  0  1^ 
'  z _ o' 


CP 


in 


(64) 


(63) 


where  c  is  the  intrinsic  impedance  of  Che  medium, 
power  to  Che  array  is 


P 


in 


1 

2 


I 


The  time-average  input 

(66) 


Since  nonzero  voltages  exist  only  at  the  feed  points,  (66)  can  be  written 
as 
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(67) 


in 


i »» (i")!) 


i  Ite  (v‘y*V) 

* 

Y  +  Y 

1  +  il  '*1 

±  V  (-= - V 

2  ~1  ^  2  ^ 


where  is  a  column  matrix  reduced  from  V  by  retaining  only  the  N  (nuo^er 
of  wire  antennas  in  the  array)  nonzero  feed-point  voltages,  and  Y^  is  an  N«N 
admittance  matrix  reduced  from  Y  by  deleting  all  elements  that  do  not  cor¬ 
respond  to  the  feed  subsections. 

In  the  far  zone,  the  magnetic  potential  can  be  approximated  by 

M 

^z  "  exp(-Jkr)  I  I  exp  I  Jkr  cos  a  )  (68) 

q-1  H  s  4 

where  is  the  angle  between  the  position  vector  to  the  qth  subsection 
carrying  current  and  that  to  the  field  point.  Using  the  matrix  repre¬ 
sentations  of  F  as  in  (6)  and  of  I  as  in  (62),  we  can  write  (68)  as 

\  ««P<-Jkr)  F  I  (69) 

or  in  reduced  matrices  as 


For  simplicity  we  write 


5i 


Y 


With  (67) g  (70)  and  (71),  D  In  (65)  becomes 


(ojliA£  sin  0  )^  v!  G,  g!  V, 
o  -1  ~1  ~1  ~1 


Vi  I 


Y  +  Y 


)y, 


(71) 


(72) 
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It  can  be  shown  [23]  chat  Che  excitation  voltage  matrix  required  for  making 
D  In  (72)  a  maximum  Is 


and 


Xi  ♦  -1 

!!i  •  Si 


2  Y  +  Y 

_  (uiuAQ.^  f  rk_iiic 
°M  4rn  £l  I  2 


(73) 

(74) 


The  directivity  optimization  problem  with  voltage  excitation  is  now  wom- 
plcccly  solved,  and  the  effect  of  mutual  coupling  has  been  taken  into 
consideration.  As  can  be  seen,  the  main  cask  in  obtaining  numerical 

solutions  lies  in  the  determination  of  the  admittance  matrix  Y, . 

~1 

Using  Che  above  procedure,  Cuianins  [23]  determined  the  maximum  di¬ 
rectivity  in  Che  princ^.pal  H-plane  of  a  circular  array  of  4  uniformly  spaced 
center- fed  wire  antennas.  Figure  7  shows  the  variation  of  maximum  directivity 
(0^)  versus  antenna  length  (2h/X)  with  array  diameter  (d/A)  as  the  parameter. 
The  points  marked  by  crosses  correspond  Co  Che  maximum  directivity  of  a 
circular  array  of  4  isotropic  sources  as  computed  by  the  method  outlined 
in  Section  III. 

The  solution  of  Che  dual  problem  of  an  array  excited  by  a  set  of 
current  sources  follows  an  entirely  similar  procedure  [23]. 
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VIll.  OTHER  CONSIDERATIONS 


We  discuss  here  several  related  aspects  of  the  array  optimization 
problem  which  either  can  be  treated  by  an  extension  of  some  of  the  pre¬ 
ceding  techniques  or  need  special  attention. 

(a)  Maximization  of  Power  Gain  -  When  the  ant'snna  wires  are  not 
perfectly  conducting,  resistive  losses  occur  and  power  gain  is  no  longer 
the  same  as  directivity  or  directive  gain.  Electromagnetically  speaking, 
in  (61)  on  the  surface  of  subsection  p  is  no  longer  zero,  but  is  equal 
to  the  product  of  1^  and  Z^,  the  internal  impedance  per  unit  length  of  the 
wire  conductor. 


where 


and 


V  -  E  (A£)  -  Z*  I  , 
P  *P  PP  P 


Z'  -  Z,(At) 

PP  1 

Z  . 

i  '  0  ^  o6 


(75) 

(76) 

(7/) 


In  (77),  0  and  6  are  respectively  the  conductivity  and  the  skin  depth  of  the 
wire  conductor  at  the  operating  frequency.  Hence,  if  the  generalized  im> 
pedance  matrix  Z  for  perfectly  conducting  wires  has  been  found,  the  only 
modification  needed  for  a  moment  solution  with  finitely  conducting  wires 
is  the  addition  of  Z'  ■  Z. (Ai)  to  the  diagonal  elements  of  Z. 

PP  1 

It  is  obvious  that  the  achievable  maximum  power  gain  for  a  given  array 
configuration  is  lower  than  the  maximum  directivity  because  of  the  resistive 
losses.  Numerical  results  have  also  shown  [23]  that  for  arrays  with  closely 
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spaced  elements  the  required  excitations  for  maximum  power  gain  are  not 
the  same  as  those  for  maximum  directivity. 

(b)  Arrays  with  Random  Errors  -  Under  practical  applications  random 
errors  exist  in  excitation  amplitudes  and  phases  as  well  as  In  element  posi¬ 
tions.  It  is  then  of  Interest  to  examine  the  effect  of  random  errors  In 
these  design  parameters  on  the  optimization  procedure.  This  has  been  done 
for  arrays  with  an  arbitrary  geometrical  configuration  [27-30].  Correla¬ 
tions  are  allowed  to  exist  between  the  errors  In  the  array  parameters  and 

no  restrictions  are  necessary  either  on  the  magnitude  or  on  the  probability 
distribution  of  the  random  errors.  The  dependence  of  the  expected  direc¬ 
tivity  or  SNR,  the  main-beam  radiation  efficiency,  the  optimum  excitation 
amplitudes  and  phases  and  the  radiation  pattern  on  the  variance  and  cor¬ 
relation  distance  of  parameter  errors  has  been  studied  [29,30].  It  was 
found  that  the  excitations  calculated  on  the  basis  of  no  random  errors  do 
not  yield  a  maximum  expected  directivity  when  parameter  errors  exist. 

(c)  Techniques  for  Large  Arrays  -  In  the  numerical  solution  of  the 
array  optimization  problem  by  the  method  of  moments  In  Section  VII,  It  Is 
necessary  to  invert  the  generalized  impedance  matrix  Z,  as  defined  in  (62). 
The  order  of  the  matrix  Z  is  M«NS,  which  Is  the  product  of  the  number  of 
wire  antennas  In  the  array  and  the  number  of  subsections  for  each  antenna 
(assuming  equal  antenna  lengths).  Hence  M  can  be  very  large  for  large 
arrays  with  many  elements,  especially  when  the  elements  are  of  a  length 
which  is  an  apprecl.tLle  fraction  of  the  operating  wavelength.  The  feasi¬ 
bility  of  inverting  large  matrices  Is  constrained  by  computer  memory 
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capacity  and  the  cost.  Alternative  techniques  which  will  relax  these  con¬ 
straints  are  therefore  of  great  importance. 

For  circular  arrays  with  many  uniformly  spaced  elements,  the  Inver¬ 
sion  of  the  relevant  large  matrix  Is  facilitated  through  dlagonallzatlon 
by  a  change  of  basis  with  a  unitary  transformation  matrix  [31],  The  columns 
of  the  transformation  matrix  are  th.t  eigenvectors  of  a  rotational  operator. 
With  this  technique  the  Inversion  of  the  large  matrix  can  be  evaluated  by 
straightforward  matrix  multiplication. 

For  array  configurations  with  no  rotational  symmetry  other  methods 
must  be  sought.  The  far-zone  electric  field  due  to  an  array  of  N  parallel 
z-dlrected  dipoles  Is 


E(e.<^) 


-J6  r  N  j6  r  ‘Q 


N  IB  r  -u  f 

-h 


JB  z'cos  6 

)e  sin  6  dz*  ,  (78) 


where  6  Is  the  phase  constant,  r^  is  the  vector  from  the  origin  to  the 
center  of  the  nth  dipole,  and  d  is  the  unit  vector  from  the  origin  to 
the  observation  point.  The  first  step  in  the  optimization  problem  Is  the 
Insertion  of  (78)  In  the  expression  for  the  performance  Index  of  Interest. 

In  order  to  maximize  the  performance  index.  It  Is  necessary  to  use  an  appro¬ 
priate  expansion  for  the  currents  I^(z').  The  moment  method  discussed  In 
Section  VII  used  expansion  functions  defined  over  subsections,  which  re¬ 
sulted  In  M  by  M  matrices.  An  alternative  Is  to  use  the  three-term  theory 
for  cylindrical  antennas  developed  by  King  and  his  associates  [ 32 ] .  This 
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has  been  done,  and  it  has  been  found  Chat,  by  using  some  special  properties 
of  the  integrals  Involved,  accurate  numerical  solutions  of  array  optimiza¬ 
tion  problems  can  be  obtained  by  working  with  matrices  of  order  N  (not  3N 
as  first  suspected,  where  N  »  M/S).  Details  of  this  technique  together 
with  numerical  results  will  be  reported  in  Parc  (B). 
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APPENDIX  to  PART  (A)  -  PROOF  OF  'THEOREM  2 


P  -  o  -  2x'e  +  x'C  X  . 

(A-1) 

which 

Is  the  denominator  of  the  quantity  In  (32).  It  Is  necessary 

to 

prove 

that  If  C  Is  positive  definite  P  will  be  minimum  at 

-1 

(A-2) 

-  C  B 

and 

min  P  -  a  -  B'C”^6  . 

(A-3) 

Proof: 

If  C  Is  positive  definite.  It  Is  known  that  [15] 

or 

(A-4) 

x'Cx  >  - ir-  (x’  B)^  , 

B'C  B 

0m  tm 

(A-5) 

where 

the  equality  sign  applies  when 

X  -  x„  -  c"^6  . 

(A-2) 

Let  c 

■  b'c”^  B  >  0,  and  b  •  x’  B.  We  have,  from  (A-1)  and  (A-5) 

2 

P  ■  A-2b  +  x'Cx  >  A  -  2b  +  —  . 

(A-6) 

But 

2 

A-2b  +  ^-  A-  c+  7  (c-b)^  A  -  c  . 
c  c  ~ 

(A-7) 

Combining  (A-6)  and  (A-7),  we  obtain 

?’  A  -  B'c'^6  , 

(A-8) 

where  the  equality  sign  holds  with  (A-2);  hence  theorem  2  Is  proved. 
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PART  (B).  INTEGRAL-EQUATION  APPROACH  FOR 
OPTIMIZING  ARRAYS  WITH  MUTUAL  COUPLING 


I.  INTRODUmON 


In  Parc  (A)  we  have  discussed  various  techniques  for  optimizing  some 
chosen  performance  indices  of  antenna  arrays.  When  Che  array  elements  are 
parallel  wire  antennas,  Che  method  of  moments  can  be  used  for  optimization 
which  includes  the  effect  of  mutual  coupling.  As  explained  in  Part  (A)  - 
Section  VII,  it  was  found  convenient  to  incorporate  Che  method  of  sub¬ 
sections  Co  convert  Che  governing  integro-differential  equations  into 
matrix  equations  by  using  pulse  expansion  functions  and  impulsive  weight¬ 
ing  functions,  For  an  N-elemenC  array  each  subdivided  into  S  segments,  an  M  by 
M  (M"NxS)  generalized  impedance  matrix  results  which  must  be  inverted  ,..n  the 
optimization  procedure.  This  inversion  process  presents  practical  diffi¬ 
culties  when  N  and  S  are  large,  because  of  limitations  in  computer  memory 
capacity  and  in  allowable  cost.  In  this  Part  (B),  we  present  a  new  approacii 
for  array  optimization  with  the  consideration  of  mutual  coupling  that  re¬ 
quires  Che  inversion  of  only  an  N  by  N  matrix.  Since  computer  time  (cost) 

required  is  proportional  Co  Che  cubic  power  of  Che  rank  of  Che  matrix,  we 

3 

achieve  a  saving  by  a  factor  of  S  .  If  S  equals  10,  for  example,  this 
means  a  1000-fold  saving,  a  considerable  factor  indeed. 

The  approach  we  Cake  starts  with  an  integral-equation  formulation  for 
Che  array  in  terms  of  Che  unknown  current  distributions  on  Che  array  elements. 
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Instead  of  using  the  methods  of  moments  and  subsections,  we  expand  the 
current  distribution  functions  as  superpositions  of  suitable  sinusoidal 
functions.  In  particular,  we  make  use  the  three-term  theory  developed 
by  King  and  his  associates  [32].  The  subsequent  theoretical  development 
Is  quit?  Involved,  but  we  have  succeeded  In  a  considerable  reduction  In 
the  order  (and  rank)  of  the  matrices  Involved.  Numerical  solutions  for 
typical  arrays  can  be  obtained  with  a  few  dollars'  worth  of  computer  time. 
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II.  INTEGRAL-EQUATION  FORMULATION 


We  consider  an  array  of  N  parallel.  z-direcCed,  center-fed  dipole 
antennas  each  of  radius  »  and  half-length  h.  The  typical  nth  dipole  is 
center-driven  by  a  delta-function  generator  of  strength  V^.  The  dipole 
•  conductor  is  aseumed  to  be  perfectly  conducting  and  Ba  <<  1.  where  6  is  the 
phase  constant.  The  integral  equation  for  the  nth  dipole  in  terms  of  the 
currents  in  all  the  array  elements  is  (32} 


N 

I 

m"i 


m  mn 


-h 


where 


30  ‘  n 


-JBr 

K  (2.2')  -  ^ 

mn  r 


mn 


..x2 


mn 


on 


mn 


V 

:  cos  Bz  +  sin  B | z|  ] 

n  i 

(79) 

mn 

i 

(80) 

•  b  1I/2 

mn 

(81) 

m-n 

(82) 

mi*n . 

Letting  z*h  in  (79),  we  have 

N  f**  4  V 

I  I  1  (z')K  (h,z')dz'  •  -  ^  IC  cos  Sh  +  ^  sin  6h) 

^ ,  j  n  osi  ju  n  £ 

in*l  i 
-n 


(83) 


Eliminating 

N 

I 

m*i 


C^cos  6h  from  (79)  and  (83),  we  get 
h 

r 

I  I  (z')[K  (z,z')  cos  8h  -  K.  (h,z')  cos  Bzjdz* 

m  mn  mn 

-h 

•  60  '^n  ® 


(84) 
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Equation  (84)  can  be  rewritten  as 


I  (*')K'  (z.z*) 
n  mn 


cos  8h  dz' 


where 


and 


■  ^  (V  sin  6  (h-|z|)  +  u  (cos  8z  -  cos  6h)]  , 
ou  n  n 


K'  (z.z*)  -  K  (z,z')  -  K  (h,z') 
mn  *  mn  *  mn  * 


N 

-J60  I  I  (z')K  (h,z')dz' 

m  mn 

m«l 


With  n«l,2,...,N,  we  obtain  N  simultaneous  Integral  equations  from 
(85). 


(85) 

(86) 

(87) 
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ill.  THE  THREE-TERM  THEORY 


King's  chree-cerm  theory  approximates  the  currents  Ig(z)  with  three 
parts:  one  part  is  a  sinusoid  maintained  directly  by  the  driving  voltage, 
and  the  other  two  parts  are  a  shifted  cosine  and  a  shifted  cosine  with 
half-angle  arguments  which  are  Induced  by  coupling  between  different  parts 
of  the  antennas. 


where 


tD 


Sj^(2)  •  sin  S  (h  -  |z| ) 


(88) 

(89) 


82(2)  ■  cos  6z  -  cos  6h  (90) 

Sj(z)  ■  tos  y  6z  -  cos  j  6h  (91) 

The  Integral  on  the  left  side  of  (85)  possesses  the  following  approximate 
properties  for  the  currents  In  (88)  -  (91)  for  different  ranges  of  Bb^ 
values  (33): 

(a)  For  5b  1, 

mn 

h 

[  I^(z')K^j^(z,z')  cos  Bh  dz'  * 

-h 

where  the  function  K'  denotes  the  real  part  of  K'  (z,z'). 

mnK  nn 

(b)  For  3b^  j.  1, 

mn 

h 

f 

I  Sj^(z')K^j^(z,z')  cos  Bh  dz'  3  1^2  §2(2)  ■  (93) 

-h 
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(c)  For  all  values  of  Bb 


mn* 


S2(z')KjJjj^j(2,z' )  cos  Bh  dz'  -  i>2S2(z) 


S3(z')K^r(z,z')  cos  Sh  dz'  s 


(94) 


(95) 


h 

1  (z')K'  -(z,z')  cos  Bh  dz'  ■  C'!S,(z),  (96) 

n  nni  j  j 

-h 

where  the  function  K'  ,(z,z')  denotes  the  Inaglnary  part  of  K'  (z,z'). 

onl  mn 

The  proportionality  constants  and  in  approximations 

(92)  •  (96)  are  best  determined  where  the  distribution  functions  In  the 

integrands  are  at  their  maximum  values. 

Under  normal  circumstances  Bb _  >  1  for  mt(n  and  Bb  -  Sa  1. 

mn  nn 

SiJ)stltutlng  (88)  -  (91)  in  the  left  side  of  (85)  and  using  (92)  -  (96), 
we  have 

(a)  m-n,  •  Ba  <  1. 


h 

f  M) 


-h 


A'  ’  S.(z')K'  (z,z')  cos  Bh  dz' 
n  1  nn  ’ 


a(1)  ^i(l)  s  (z)  +  S.(z) 

n  nnR  1  n  nnl  3 


(97) 


-h 


A^^^  S-(z')  K'  (z,z')  cos  Bh  dz' 
n  i  nn 


=3<^) 


(98) 
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f  A'"'  S-(z')K'  (z.z')  cos  6h  dz' 
j  n  j  nn 


(3) 


0'^'^  S,(z) 
n  nn  3  • 


(99) 


(b)  mjtn,  db  ■  6d  >  1. 

jm  nm  — 


h 

^  S, (z’)K'  (z,z')  cos  6h  dz' 

ID  1  OUl 


-h 


h 

I  A^^^  S-(z')K'  (z,z')  cos  Sh  dz' 
in  z  mn 


j 

-h 


•  *2<'»  -  *'ll 


(101) 


h 

A^^^  S,(z')K'  (z,z')  cos  6h  dz'  -  A^^^  ij«' ^^^S,(z). 
m  3  mn  ’  m  mn  3 

-h 

The  proportionality  constants  in  (97)  -  (102)  are: 


(102) 


r*’  cos  liR,  cos  BR» 


,(U 

'mnR 


-h 


cos  6h 


II 

[  S^(z')( 


cos  BR.  cos  BR» 

— r - i - - - ^]dz'  ,  Bh  >  n/2 


-h 


(103) 


with 
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(104) 


mnR 


^.(1) 


,.(2) 

''onl 


..(3) 


*^1 

•^3 

cos 

gh_ 

1-C08 

Bh 

COS 

gh_ 

1-COS 

Bh 

-cos 

ph 

1-cos 

Bh/2 

-cos 

ph 

1-cos 

Bh/2 

cos 

Bh 

1-cos 

6h/2 

(z-2  *  1,2 

isn 

l(h-z’)^  + 

mn 

((h-2'-X/4)^  + 

nn 


cos  6R.  cos  6R« 


-h 


h  ^ 


cos  6R^  cos  BR2 


Jdz' 


sin  BR.  sin  BR, 

-  -R— ^1“^’ 


-h 


f  sin  BR,  sin  BR, 

J  V^’>l-R--^-  — 

1. 


-)dz’ 


f  ■'"•'I 

1 

-h  ^ 


-JBR.  -JBR- 


')dz' 


(105) 

(106) 


(107) 


(108) 


(109) 


(110) 


(111) 


When  6h  ■  n/2,  expressions  (109)  -  (111)  become  indeterminate  and  com¬ 
ponent  current  functions  different  from  those  given  in  (89)  -  (91)  must 
be  chosen.  This  special  case  will  be  discussed  in  Section  V. 

We  an*  now  ready  to  substitute  the  above  in  (85)  and  equate  the 
coefficients  of  S|^(z),  kal,2,3.  We  obtain  3  sets  of  equations  (nBl,2, . . .  ,N) 
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(112) 


^  m  "^ninR  60  ''n 


m«l 


oftn 


60 


m«l 


In  (113), 


“  I 

-  j60  I 

tn"l  / 


I  (z')K  (h,z')dz' 

m  ran  ’ 


npl 


m  mn 


m  mn 


where 


l-^^h) 


II 

f  S^(z')K^(h,z’)dz' 


-h 


S,(z’)K  (h,z’)dz' 

z  mn 


-h 


II 

I  S2(z')Kj|^(h,z*)dz' 


With  (115)  -  (118),  we  can  rewrite  (113)  as 


(113) 


(114) 


(115) 


(116) 


(117) 


(118) 
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(119) 


N 

I 

ni*l 


A^^^Kl  - 


m-l 


-  'I'on^h)) 


N 

I 

m^l 


where  6  is  a  Kronecker  delta.  For  an  N-element  array,  n-l,2 . N, 

on 

and  each  of  (112),  (119),  and  (llA)  represents  N  simultaneous  equations. 
It  Is  therefore  convenient  to  use  a  matrix  notation. 
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IV.  MATRIX  EQUATIONS 


Defining  N*!  column  nuiCrlces  [A^^^],  and  [V],  we 

can  write  (112),  (119),  and  (114)  as 


(120) 


-  0.  (122) 


Equations  (120)  -  (122)  can  be  Inverted  and  rewritten  as 


-  (P^^hlVl 

[A<2)]  .  (p<2)j(v] 

(a(3)]  .  (p(3))(v]  , 


whe  re 


J-  r„'^2),-i 

60  ^^mnR  ^ 


(P<2)i  .  j^(b)j-lj^(a)jjpvl)j 


and 


ip”>i 


u‘^'i 


(123) 

(124) 

(125) 

(126) 


(127) 

(128) 

(129) 
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The  current  distribution  in  the  mth  dipole  in  (88)  can  also  be 
extended  to  the  N-eleoent  array  and  written  in  a  matrix  form. 

[I]  -  S^(z)(A^^^l  +  S2(2)[A^^^]  +  S3(z)[A^^^]. 

Combining  (123)  -  (125)  with  (130),  we  obtain 


[II  -  [VIIV] 


where  the  nnth  element  of  the  admittance  matrix  [Y]  is 

+  S,(z)pJ^^  +  S  (z)p2^  . 
inn  1  Dn  £  nil  j  nm 

With  (131),  we  are  finally  ready  to  attack  the  optimization  problem. 

But,  before  we  do  this,  we  shall  take  care  of  the  case  for  6h  -  it/2 

(half-wave  dipoles)  which  would  make  (109)  -  (111)  indeterminate. 


(130) 


(131) 


(132) 
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V.  HALF-WAVE  DIPOLES 


It  is  obvious  that,  whsn  Bh  ■  ii/2,  we  could  not  have  proceeded  from 
(85)  as  we  did  in  Section  III.  We  must  start  from  the  basic  Integral 
equation  (79).  When  z  Is  set  to  equal  h,  (83)  becomes 

•h 

Subtracting  (133)  from  (79),  we  get 


N 

I 

m"l 


I  (z')K'  (z,z')dz' 
m  mn  * 


.  V 

-  ^  6z  -  (sin  6| 


z|  -  1)],  (134) 


-h 


where  K*  (z,z')  has  been  defined  In  (86).  C  can  be  found  from  (134) 
mn  n 

by  setting  z  -  0  in  (79). 


C 


n 


I  (z’)K  (0,z’)dz’ 

la  inn 


Combining  (134)  and  (135),  we  obtain 


(135) 


I  (z')K'n(z.z')dz' 
m  mn 


.  J- 


60  -  I)  +  cos  Bz)  , 


where 


n*l,2, . . . ,N 


U' 

n 


N 

J60  I 
m"l 


-n 


I  (z')K  (0,z')dz'. 

m  mn 


(136) 


(137) 


Equations  (136),  Instead  of  (85),  must  now  be  solved. 
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Thr  three-tem  expansion  for  current  distribution  in  (88)  -  (91)  must 


also  be  modified.  We  vrite 


ID 


k-l 


S'(z) 
m  k 


with 


(138) 


8^(2) 


•  sin  S  z  -  1 


(139) 


S^U)  ■  cos  Sz 


(140) 


S2(z)  •  cos  j  6z  -  cos  ^ 


(141) 


Following  the  same  procedure  outlined  in  Section  III,  we  obtain,  instead 
of  (112)  -  114), 


?  b“>  v'<1>  .  -.l-v 

m  ^mnR  60  n 


(142) 


In  (143) 


na^n 

m*l 


..(1) 

^mnR 


f  cos  6R^  cos  BRn 

-  J  4<x'>1-T— -  — 

-h  ^ 


•Idz' 


(145) 


.(2) 

'mnR 


n 

I  n 


cos  BR,  cos  BR_ 

(z')l— S - - -Uz’ 


-h 


(146) 
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,.(3) 

'mnR 


-h  '  h 


C147) 


x'  -  —  -1 

mnl  1  _  gQg 


774  / 

-h  1  ^9 


rT  Jdz'  (148) 


X-^2)  _ -1 

^ranl  - 


cos  7r/4 


6Ro 


-h 


R 

2 


.(3) 


mn 


i  “  cos  7r/4  j  ®3(z')[— 


-jSRi  -jgR 


where  R^  and  R^  have  been  defined 

h 


-h 


■]dz 


C150) 


in  (104)  and  (105) .  Also, 


where 


>^m«» 


-h 


(2)  r 

^mn  ■  j  S2^z')K^jj(0,z')dz' 

“h 

(3)  ^ 

-  /  S’(z')K  (0,z')dz' 


-h 


(152) 


(1"3) 


(154) 


In  matrix  form,  (142)  -  (144) 


can  be  solved  to  give 
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(155) 

=  [Q^^^][V] 

(156) 

=  [Q^^^][V]  , 

(157) 

where 

]  =  _  i_  (1) 

^  60  ^^mnR  ^ 

(158) 

(159) 

[Q<”l  - 

(160) 

and 

- 

[x^^^(0)][x'^^^]~^[x'^?^]  +  [x’^o^  -  X^^^(O)] 

^^mn  ^^mn  ‘  ^mnl  ^  ^*mnR  *mn 

(161) 

• 

-  [(1  -  A  )x'^D^  -  X^^^(O)]  -  [x^^^(0)][x’^^^]”^ 
mn''^mnR  ^mn  ^  '■^mn  ^  '^■^''^mn  ‘ 

lx'  ] 

^^mnl  ^ 

Combining 

(1A9)  - 

(151)  with  (138),  we  have 

(162) 

[I]  =  [Y'llV]  , 

(163) 

where 

(z) 

mn 

S,'(z)Q^^^  +  Sl(z)Q^^^  +  S;(z)Q^^^ 

1  ^mn  2  '^mn  3  ^mn 

(16A) 

Equations 

(163)  and  (16A)  for  half-wave  dipoles  correspond  to  (131) 

and 

(132)  respectively. 


-53- 


VI.  MAXIMIZATION  OF  DIRECTIVITY 


The  far-zone  electric  field  of  an  array  of  N  parallel,  z-directed, 
center-fed  dipoles  is 


n=l 


jer 


n 


T  t  i\  j3z'cos  0  .  Q  j  , 
I  (z'^e-^  sin  6  dz' 


(165) 


-h 


where  S  is  the  phase  constant,  r^  is  the  vector  from  the  origin  to  the 
center  of  the  nth  dipole,  and  0  is  the  unit  vector  from  the  origin  to  the 
observation  point.  Equation  (165)  can  be  rearranged  and  written  in  a 
matrix  notation. 


E(e,(^) 


[H]^I]  dz'  , 


(166) 


where  both  [H]  and  [I]  are  N^l  column  matrices. 

sin  6  _  ^  ^n 


[H]  is  H  : 
n 


H 


4Tr 


(e 


The  typical  elements  of 

-jSz'coBe,  ^ 


Substituting  (131)  in  (166),  we  obtain 


where 


E(0,<J.) 


h 

[Hj’^lYjdz' 

. 

-h 


[Mj'^V]  , 


[M]' 


h 

* 

[H]^[Y]  dz' 
-h 


[VJ 


(168) 

(169) 


54- 


is  a 


If  Sh 


where 


IxN  row  matrix.  Of  course,  (163),  instead  of  (131),  would  be  used 
=  it/2.  The  only  effect  would  be  to  change  [Y]  in  (169)  to  [Y']. 
Now,  the  array  directivity  is,  from  (65), 


r^|E(e  ♦  )|2 

-  --y-  — 

in 


(170) 


is  the  time-average  input  power  to  the  array. 


=  i  Re  {[V]  lYj^^^Q-[V]} 

=  7  LVl'^[Y^][V].  (171) 


The  elements  of  [Y^]  are  the  real  part  of  those  of  the  driving  point 
admittance  matrix  (168)  and  (171),  0(9^, (jj^)  in  (170)  becomes 

D(9  ,^  )  - - 2 - 2 - 

°  °  30  IV]^[Y^]  [V] 

where  [M^]  is  (M]  in  (169)  evaluated  for  the  direction  (6q»'}’q)»  Equation 
(172),  as  (7),  is  a  ratio  of  Hermitian  forms,  and  Theorem  1  in  Section  III, 
Part  (A),  can  be  used  to  find  the  maximum  directivity 
required  voltage  excitations  [Vj.  We  have 

“m-  M 
and 
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The  optimization  problem  is  now  completely  solved.  The  required 
matrix  [Y^^]  is  obtained  easily  from  the  admittance  matrix  [Y]  (or  [Y'] 
if  Sh  =  "72),  whose  formulation  has  been  developed  in  the  previous  sec¬ 
tions.  We  note  that  [Y  ]  is  of  a  dimension  NxN  for  an  N-element  array 
irrespective  of  the  length  of  the  dipoles. 
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VII.  NUMERICAL  EXAMPLE 


The  integral-equation  approach  formulated  above  for  optimizing  arrays 
with  mutual  coupling  is  applied  to  a  four-element  circular  array  for  which 
some  important  results  on  directivity  optimization  have  been  obtained  by 
using  the  method  of  moments.  Although  this  approach  applies  to  larger  arrays, 
it  was  thought  advisable  to  check  with  some  known  results  first  in  view  of  the 
rather  involved  analytical  process.  Cnee  the  formulation  and  the  results  have 
been  verified  for  the  four-element  array,  the  extension  to  larger  arrays  is 
straightforward  and  needs  only  slight  changes  in  the  computer  program. 

The  array  has  four  parallel,  z-directed,  dipoles,  each  of  diameter  2a 
and  length  2h,  uniformly  spaced  around  a  circle  of  diameter  d.  The  coordinate 
system  is  selected  such  that  one  dipole  coincides  with  the  z-axis,  a  second 
one  lying  in  the  xz-plane,  and  a  third  one  lying  in  the  yz-plane.  The  follow¬ 
ing  parameters  are  chosen: 

2a/A  »  0.0025 
2h/A  =  0.36  (h/a  «  144) 
d/X  *  0.61  . 

The  directivity  of  tius  array  in  the  principal  H-plane  (0  =  90®)  is 
maximized  for  each  value  of  (ji  by  adjusting  the  amplitudes  and  the  phases  of 
the  excitation  voltages  V^,  n<*l,2,3,4.  The  results  are  plotted  as  the  curve 
in  Fig.  8,  which  coincides  almost  exactly  with  that  obtained  by  Cummins  [23], 
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who  used  the  method  of  moments  and  the  method  of  subsections  described 
In  Section  VII,  Part  (A).  The  Fortran  program  used  for  computing  the 
maximum  directivity  In  Fig.  8  Is  appended  to  show  Its  relative  simplicity. 
In  spite  of  the  complexity  of  the  formalism,  the  cost  of  computing  the 
entire  curve  In  Fig.  8  on  an  IBM  360/50  computer  was  only  about  five 
dollars . 
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VIII.  CONCLUSION 


The  obvious  advantage  of  the  Integral-equation  approach  In  solving  the 
directivity  maximization  problem  for  an  array  of  dipole  antennas,  as  compared 
to  the  method  of  moments  using  subsections,  Is  the  S-fold  reduction  In  the 
order  of  the  matrix  to  be  Inverted,  where  S  Is  the  number  of  subsections  for 

each  dipole.  This  results  In  a  saving  In  computer  time  (cost)  by  a  factor 

3  * 

of  S  ,  In  addition  to  relaxing  the  requirements  on  memory  capacity.  Besides 

the  directivity,  quantities  such  as  the  current  distribution  on  each  dipole, 

the  self  and  mutual  Impedances,  the  radiation  pattern,  etc,,  can  all  be 

calculated  without  difficulty. 

For  arrays  with  many  elements  (N  very  large),  practical  computing  dif¬ 
ficulties  will  arise  even  If  the  antennas  are  not  divided  Into  subsections. 

In  such  cases,  other  techniques  are  needed  to  simplify  the  computing  procedure. 
Because  of  the  existence  of  rotational  symmetry  In  a  uniformly  spaced  circular 
array.  It  Is  possible  to  circumvent  the  necessity  of  Inverting  any  matrix  by  the 
Introduction  of  a  rotational  operator  [17].  Hence,  the  number  of  elements  In  a 
uniform  circular  array,  no  matter  how  large,  represents  no  real  constraint  on  the 
feasibility  of  obtaining  numerical  solutions.  On  the  other  hand,  linear  arrays 
possess  no  circular  symmetry:  thus  the  technique  of  using  a  rotational  operator 
does  not  apply  and  other  methods  must  be  sought  when  N  Is  very  large.  Special 
methods  for  handling  pattern  synthesis  and  performance  optimization  of  very  large 
arrays  constitute  an  Important  area  for  further  research. 

For  the  particular  example  In  Section  VII,  symmetry  about  the  plane  (xy-plane) 
bisecting  the  dipoles  reduces  the  effective  order  of  the  matrix  handled  by  the 
method  of  moments  by  a  factor  of  2.  Symmetry  property  about  the  plane  contaln- 
Ingthe  diametrically  opposite  elements  simplifies  the  computation  for  both  the 
moment  method  and  the  Integral-equation  method. 
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APPENDIX  TO  PART  (B) 


FORTRAN  PROGRAM  FOR  COMPUTING  DIRECTIVITY  OF  A  FOUR-ELEMENT 

CIRCULAR  ARRAY  AS  A  FUNCTION  OF  AZIMUTH  ANGLE 

)  , 

S'  iMrt'  II 1 !  1  .'It  (  A1 ,  ‘^r.  f ‘<C  ,H  ,■ ! .  L  ) 

r.iii  PI  i-y =-n  '  n (i-i.i. )  .'•r, (t.  .«*• )  .^r. (.>i .n  j 

c 

’rnT” 

Dll  1<'  .1-1  ,  ' 
i>r.(  I  ..i)=(n.  ,0.) 

c 

U) 

«C(  1  .  J)=>-T.(  1  I  ,K)-KC(K.  J) 

Pf-Tl'P'i 

t;''D 

stir  ntiin  IT-  r.Ai.-'-'M 

r.!l!  PI.!  ..M'  (••Srt)  ,MI  ,rt:-.(4,4  )  ,  S  1  ( ,'i  )  «  1  1  )  t  Kri  (‘t  ) 

<. 

1- /\I.  •• 1  1  •  (. -!  t  f  .  .'i  * 

D 1  =  ( n . ,  1  .  ) 

|l(l  IHA  I  =  1.M 

V 

1  fit) 

TTnT‘’“".rnr.  '•• 

!r  (  1-.!)  « 

D  '(  I .  .1  )  =  (n.  ,o.  ) 

>  t 

h)':! 

(.•ri  1 1> iv  1  *  .1  , 

C 1  =?.'.=  1  .f’n(-n*Hlj;iV  (  (•U-l)>5(  !-l>  »/hLUAT  (i  i ) -kA  .  I**  1  s'-)? 
p.i'l  I  1  I  ,.l  )  +  '•'(  1  .it  F  I'MC'S  (Cl  (C)  )  ) 

<■ 

ir>/i 

TTjTi ,  .rj  =  .i  FT ;:n 

c.i)  n*  106 

Biv(  I  ,.i)  =  io.,n. ) 

o 

1. 1  * 

TJTiTJTTpr 

IH)  1  I-l,"' 
nn  1  .isi,'' 

c 

inv 

“C? -TtlTT- '((  '')  ■  '  ' 

C4=si>pf  ) 

CVs?.'"! ,  noon^Hi  ii.AT  ( 1-1  1  )/rLuA.f  (u) >i:^ . ) 1  :>4;;h33 

1 

l 

“STiffTi'^t  K- 1  d"  t !“.?  n  /x~fi 

11  ( j,  1  )  =  (r.".s(C2)-i’.i*.siii(c?)  )/(>> 

CiI('TI:MIi; 

.1 

L- 

-ri'-T-i.  n  “T  ri  -T^ - 

"(1  H'  .1  =  1,1'  r\9>^ 

PMC  1,.|)  =  (0.  ,n. ) 

in 

TIT ‘in 

Ri'l  I  I  f  J)+S1  (  I  .O^KMCK,.!) 

fill  1 1  t  =  1  .  iM  t''' 

c 

•lii  11*  .r='i,i'' 

A«(  I  ,.n  =  ('''.  ,0. ) 
on  11  !<=i,l-i 

c 

- IT 

"Tomrm  f  u  ,r.  i 't,- ni. ,  .rj— 

RbTlIH  ■' 

L 

-TTI  T-v  -IT.P.SV  !'■' 

Ctii'Pl.r  >'  '•«  AA  (A  )  .  •  I'  (^  ,'‘)  ,AC  (A  .  All  (<»  ,4  )  ,  .''ri  (H  )  , 

1AB(4.A).m/' (A,4)  (4.A) 

- [71*^“  — !')  I  ,•  r,  r*;  ,‘*Tr  M  » '  r  < i  »  i  i  i 

''  ),  .’M  (  ^  I  (  A  )  ,t<.l  (<•.<•)  ,f<  ■<(<•,/.).  V  .P'i  (<.  )  1 

L 

. 

TTSr'T’i  'TT TL  ,  A  1  ••  •  1  ■ . '  •  1  P  »  1  J  f  r  1  f  )  ,  1.  1  M  ,  i  I  .  f  J  ,  >  '  \  a  ,  i  »  ,  >  •  i  »  v  v  i  i  , 

1PM<\S  (•'• )  .  )  ,(.A  1  I  ,A-  S,S«l  (-•.  1  ,VT  (B  ) 

Cili'PLtyv**  YV  (4,4  )  ,04,07,  Ji)(4  )  ,'0:  (4  )  ,  YF  (4  ,<• )  ,V  1  ,V  (‘  )  |1;a  In 

U 

IT  =  ('■.,  !  .1 

MLfci:l=o.i‘< 

PsO.AI 

-TM  .1 .  1  1  =  ■ 

0(2,1  )=p-"’.707 

||(  1 .1  Fs-'i.r  T /!<•/, . 

- T;rr7  r;i  =  \ ^ 

!•=  I  I  -  I  ) 

/l-L"Ar  (ir ) 

|..  1  =  11-  \ 


SSd  ,:c  +  l)  =  (HLriAf  (' 
3  TT(  I  ,'^  +  l  )  =  (FL-'.'T(<' 


••  II.  1  t  - 1  I  1  ■-  I--  ( 

7  SR(K1.ti!?)  =  .SS(liA)* 
on  1=1.4 


i''=i  A«s(  i-.n 
IF(  n-1  )SH,S4,*^4 


4.?K3? 


P'H  1  S'MLt  I  J-m-i'SMfl.,?)  »/«»b 

i;i)  TM  f>7 


MKPOJ*  |Ai>3t*oii.7n  /c4;i.-v,fwKrin| 

n*!  7'»l  I*<  .* 

••fl3  l'».  <  t).  I 

UM  70?  .1=1,* 

UXel'  .1) 

■7iI7*TT/:  STT^  riTHn  M*  vrr^TT.VJ » ^  1^1  I'.r. !  J 


II.IISII'-. )  «  FusI  li  I**!  l.l.T)  J?,-l 

J,  I  Fv?.o(?.wnH-AniiiuM7  J  ) 

|tY*««  I  ■■  HO.  viiA 


.-■••j  I « r;*  ( j ,  i )  I 


701  •  f  I  IF-  j  H,  7or..»)  I  ....  (  I  I ,  I  .I..>n  I 
iiiiTT  {•  I a”  I**”. •  * 7** li"  - ^  ^ »  I  *  ,  *  )  =  '  •  I 
C/LL  CAI  •VMYY.y  I  ,*J 

pi  c»f.n  1  =  1.* 


I 

PfI  Hf"!  .1=1,* 
wr.l  wi=\/n.v|  (  I, .Ip-  I  -i.ij 

- rrm  =-^ - 

vvj  I  j=c/‘'<;<v) ) 

V:i=^^AL  (  VI  ) 

- vrr— 7,1^3  crm - 

V!.=v/./y^ 

P«AS(  I  )=-AlA  (V'J'Jlf  O./H.)*!* 

-- RTT-’-ff  TTrTT;7^~i^T-.-nT7T7-r~i-m"rT'',  -»Ti7;Tti , - 

PHOO  hiJUriM  (  •  O'  ,  '  V(  '  ,  I  H,  '  J  =  •  ,?.••  1  ?.^/?X,  •  VV«  •  ,  IH  ,  •  )  =  •  ,c  1  ?  .?/?A  ,  •Khfti;  (  • , 

n  3 .  •  I  -  • ,  <■- 1  ? .  •)  J 

- rniTT— roTT'VT - 

n»(  hii?  !« 1,* 

fin?  r.Aii‘=r;Ari+ci:n  )> '<(  1 1/50. 


I-/.  I  I  I ....  I  •  ) 

POH  k  1.  m  (5,  pool  JO  All 
pool  Pdfl.A  I  (  •  I  •  ,  •Cf.h'v  •  ,(•  I  ?.t?  J 


nr 


23 


ARRAY  DIAMETER.  Zpf\ 


Fig.  3  -  Directivity  and  main-beam  radiation  efficiency  for  a 
circular  array  with  isotropic  elements  (N-12,  6  -90*. 
♦o-O*)- 
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MAIN*  BEAM  RADIATION  EFFICIENCY 


T(u) 


Fig.  A  -  Spatial  distribution  of  noise  power;  u  ■  (2iTd/X)sln  6. 
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Fig.  5  >  SNR  optimization  for  7-eleinent  broadside  array. 
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DIRECTIVITY, 


ARRAY  DIAMETER,  Zp/\ 


Fig.  6  -  Optimum  directivity  and  main-beam  radiation  efficiency 
for  a  circular  array  with  dipole  elements  (N-12,  6  -90”, 
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MAIN- BEAM  RADIATION  EFFICIENCY. 


Maximum  directivity  of  a  A-element  circular  array  In 
direction  9«90*  and 
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Fig.  8  -  Maximum  directivity  of  a  4-elemcnt  circular  array  aa  a 
function  of  azimuth  angle. 
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